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I.  INTRODUCTION 


Researchers  have  used  a  variety  of  techniques  ~  to  obtain  blast  loading 
data  in  their  study  of  generic  shapes.  Of  these  basic  shapes;  rectangles, 
cylinders,  spheres,  and  cones  have  special  interest  for  the  Army  in  that 
many  weapons’  systems  include  one  or  more  of  these  shapes.  The  blast 
loading  data  were  obtained  in  the  studies  by  means  of  shadowgraphy , 
interferometry,  pressure-time  recording,  and  with  force  balances  to  define 
the  flow  fields  and  drag  loading  generated  by  the  blast  loading  from  a  given 
weapon. 

More  recently,  two  and  three  dimensional  hydrocodes  have  been  developed 
and  are  being  used^“8  to  predict  the  flow  fields  created  by  blast  wave  loads. 
In  a  sense,  the  hydrocodes  may  produce  predicted  results  that  summarize 
the  data  from  the  several  experimental  techniques  used  to  determine  blast 
loads.  Wave  systems,  flow  fields,  and  pressure  or  density  fields  may  be 
generated  by  output  from  hydrocodes.  Integration  of  the  predicted  loads 
over  a  given  target  surface  may  be  used  to  calculate  total  loads  and  to 
generate  coefficients  of  pressure  or  drag. 


1 

S.H .  J.  Allen  and  W.  G.  Vincenti,  " Wall  Interference  in  a  Two-Dimensional 
Flow  Wind  Tunnel  with  Consideration  of  the  Effect  of  Compressibility  y  " 
Nat .  Adv.  Comm .  for  Aeronautics  Report  782 y  1944 . 


9 

R.N.  Holbyer  and  R.E.  Duffy  r/The  Effect  of  Wall  Boundary  Layer  on  the 
Diffraction  of  Shock  Waves  around  Cylindrical  and  Rectangular  Obstacles y  n 
University  of  Michigan  Report  50-2 y  1950 . 


7 

W.  Bleakney  and  D.R.  White,  "Shock  Loading  of  Rectangular  Structures , " 
Dept,  of  Physics  Tech.  Report  II,  11,  Princeton  Unit.,  January  1952. 


^ N.K .  Delany  and  N.E.  Sorensen,  "Low  Speed  Drag  of  Cylinders  of  Various 
Shapes,"  Nat.  Adv.  Comm.  Aero.  Technical  Note  3038,  Wash.,  DC,  1953. 

5 George  A.  Coulter  and  William  T.  Matthews,  "Coefficients  of  Drag  Measured 
with  a  Force  Balance,"  BRL  Technical  Note  1155,  December  1954. 

Q 

M.A.  Fry  and  others,  "The  HULL  Hydro-Dynamics  Computer  Code,"  AFWL-TR-76- 
183,  U.S.  Air  Force  Weapons  Lab,  Kirtland  Air  Force  Base,  NM,  September 
1976. 


n 

Richard  E .  Lotteroy  nA  Detailed  Comparison  of  3-D  Hydrocodes  Computations 
for  Shock  Diffraction  Loading  on  an  S-280  Electrical  Equipment  Shelter y  n 
BRL  Technical  Report  ARBRL-TR-02334y  June  1981 . 

8 

John  D.  Wortman  ,  "Blast  Computations  over  a  Hemicylindrical  Aircraft 
Shelter,"  BRL  Memorandum  Report  ARBRL-MR-03115,  July  1981. 
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To  build  up  confidence  in  a  given  code  or  analytical  method,  point-by¬ 
point  experimental  comparisons  are  needed.  The  purpose  of  the  tests  reported 
here  is  to  provide  blast  loading  data  for  a  two-dimensional  cylinder  which 
will  allow  such  a  comparison.  The  loading  function  used  is  a  decaying  wave 
provided  by  the  Ballistic  Research  Laboratory's  (BRL)  57.5  cm  shock  tube. 

The  next  section  describes  the  test  procedure  used. 


II.  TEST  PROCEDURE 

The  test  procedure  is  discussed  in  three  parts:  the  BRL  57.5  cm  shock 
tube,  the  cylinder,  and  the  instrumentation. 

A.  BRL  57.5  cm  Shock  Tube 

g 

The  BRL  57.5  cm  shock  tube  was  modified  by  installing  a  92  cm  driver 
section.  Other  parameters  were  left  the  same.  The  short  driver  allowed  the 
rarefaction  wave  from  the  driver  section  to  catch  up  with  the  shock  frontlO 
and  create  an  exponentially  decaying  shock  wave  in  the  shock  tube  test 
section.  Figure  1.  The  pressure-time  history,  although  derived  one- 
dimensionally  in  the  shock  tube,  simulated  quite  well  a  free-field  wave. 
Comparisons  between  the  two  types  of  waves  are  made  in  the  Analysis  Section. 

The  shock  tube  was  operated  in  the  air-air  mode  with  ambient  air  present 
in  the  test  section.  The  entire  downstream  test  section  was  open  to  the 
outside  air.  No  end  plate  was  used.  Driver  pressures  were  chosen  to  give 
side-on  shock  overpressures  in  the  range  42.3  kPa  (6.1  psi)  to  112.2  kPa 
(16.3  psi) .  Shots  were  repeated  until  pressure-time  histories  were  recorded 
for  15°  increments  of  cylinder  rotation  for  the  pressure  range  chosen. 

B.  Two  Dimensional  Cylinder 

Figure  2  shows  the  location  of  the  transducers  and  Figure  3  shows  the 
cylinder  in  50.8  cm  square  test  section  of  the  shock  tube.  Because  of 
transducer-connector  length,  four  transducer  locations  were  chosen  at  45° 
intervals  around  the  model  near  the  center  of  its  length.  The  remaining  15° 
intervals  were  exposed  to  succeeding  repetitive  shots  of  the  shock  tube. 

Station  0  was  defined  to  be  facing  directly  into  the  on-coming  shock  front,  0°. 


George  A.  Coulter  and.  Brian  P.  Bertrand ,  "BRL  Shock  Tube  Facility  for  the 
Simulation  of  Air  Blast  Effects ,  "  BRL  Memo  Report  No.  1685,  August  1965. 

10 

C.W.  Lampson ,  "Resume  of  the  Theory  of  Plane  Shock  and  Adiabatic  Waves  with 
Applications  to  the  Theory  of  the  Shock  Tube, "  BRL  Tech.  Note  No.  1Z9, 

March  1950. 
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Figure  2.  Initial  location  of  transducers. 
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Figure  3.  Two-d i mens i onal  cylinder  installed 


The  cylinder  was  fastened  to  the  round  access  port  in  the  top  of  the  test 
section.  The  top  of  the  test  section  was  graduated  in  degrees  so  the  port 
could  be  rotated  and  matched  in  rotation  (with  the  cylinder)  for  the 
succeeding  transducer  locations.  Appendix  A  gives  details  of  the  test 
cylinder. 

C.  Instrumentation 


The  schematic  of  Figure  4  summarizes  the  data  acquisition-reduction 
system.  The  test  cylinder  contained  four  PCB  13A24,  quartz  element,  pressure 
transducers  flush-mounted  with  surface  of  the  cylinder.  Signal  conditioners 
and  data  amplifiers  transferred  the  pressure-time  histories  to  an  FM  7600 
Honeywell  recorder. 

Quick-look  data  could  be  obtained  by  means  of  the  Honeywell  1858  CRT 
Visicorder.  Final  data  processing  was  accomplished  with  a  Textronix  4051 
computer  and  related  accessories.  Final  report-ready  copies  of  pressure-time 
histories  with  engineering  units  were  made  with  the  data  system. 


III.  RESULTS 

Figure  5  shows  typical  side-on  and  stagnation  pressure  records  for 
input  pressures  of  42.3,  75.9,and  112.2  kPa. 


The  records  were  obtained  from  locations  about  0.8  m  upstream  of  the 
test  section  (Figure  1)  on  opposite  sides  of  the  shock  tube.  Ten  milliseconds 
only  of  the  records  are  shown  since  this  was  the  test  time  of  interest-through 
the  diffraction  phase  until  flow  had  been  established  about  the  cylinder. 

The  small  peaks  occurring  after  about  4  milliseconds  are  shock  reflections 
from  the  cylinder  and  test  section  walls,  and  should  be  ignored.  Table  1  is 
the  shot  log  for  the  test. 

A.  Results  for  Stations  0-30 


Pressure-time  histories  for  the  period  of  interest  are  shown  for 
comparison  in  Figure  6-18  as  a  function  of  transducer  location  and  input 
side-on  shock  overpressure.  Tables  1  and  2  summarize  pertinent  shot 
parameters  and  test  results. 

The  records  are  grouped  according  to  increasing  input  overpressure  for 
each  station  (transducer  location).  For  example.  Station  0  is  for  0°,  15  is 
for  15°,  and  so  on.  The  records  from  Stations  0-30  all  show  initially  a  peak 
of  reflected  pressure  which  is  followed  by  rarefactions  which  decay  the 
reflected  pressure  to  the  value  for  stagnation  pressure.  The  small  pressure 
peak  at  about  1.5  milliseconds  is  a  reflection  from  the  cylinder  to  the  test 
section  wall  and  back  to  the  transducer.  There  reflections  should  be  ignored. 
Similar  reflections  occur  at  later  times,  but  are  mixed  in  with  the  general 
pressure  ocillations  and  are  hard  to  distinguish. 
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15 


Figure  4.  Schematic  of  data  acquisition-reduction  system. 


PRESSURE,  kPa  PRESSURE,  kPa 


TEST:  2-0  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  5.  Pressure-time  records  for  input  waves. 
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PRESSURE,  kPa  PRESSURE,  kPa 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  5.  Pressure-time  records  for  input  waves.  (cont'd) 
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PRESSURE,  kPa  PRESSURE,  kPa 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  5.  Pressure-time  records  for  input  waves.  (cont'd) 
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28.7 

251.0 

24-81-89 

150 

23.6 

28.1 

265.7 

24-81-89 

side-on 

42.6 

28.4 

305.9 

24-81-89 

stag 

92.0/46.0 

27.2 

332.6 

24-81-90 

15 

87.7 

31.5 

330.3 

24-81-90 

60 

68.8 

31.3 

243.8 

24-81-90 

105 

41.7 

33.7 

246.4 

24-81-90 

150 

23.5 

28.9 

259.2 

24-81-90 

side-on 

42.0 

32.1 

303.1 

24-81-90 

stag 

95.6/47.0 

31.0 

335.3 

24-81-91 

30 

88.8 

30.6 

302.7 

24-81-91 

75 

61.0 

26.2 

225.6 

24-81-91 

120 

33.8 

29.2 

260.7 

24-81-91 

165 

20.0 

25.9 

255.2 

24-81-91 

side-on 

42.6 

28.5 

302.2 

24-81-91 

stag 

95.0/46.3 

27.1 

332.4 

24-81-92 

45 

80.0 

31.6 

267.8 

24-81-92 

90 

50.1 

27.7 

219.3 

24-81-92 

135 

27.0 

33.2 

265.4 

24-81-92 

180 

42.6 

28.8 

270.9 

24-81-92 

side-on 

42.0 

31.9 

298.1 

24-81-92 

stag 

94.7/45.5 

30.5 

327.9 

24-81-93 

0 

188.6 

40.0 

682.9 

24-81-93 

45 

140.6 

35.4 

509.6 

24-81-93 

90 

86.8 

43.2 

378.7 

-10.0 

24-81-93 

135 

49.0 

34. 1 

439.8 

24-81-93 

side-on 

75.5 

42.1 

582.3 

24-81-93 

stag 

174.6/92.0 

39. 1 

675.8 

*Positive  impulse  given  for  10  milliseconds. 
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TABLE  2.  TEST  RESULTS  (cont'd) 


Shot 

Transducer 

Initial 

Positive 

Positive 

Negative 

Number 

Position 

Overpressure 

Duration 

Impulse* 

Pressure 

kPa 

ms 

kPa-ms 

kPa 

24-81-94 

15 

185.4 

42.1 

665.1 

24-81-94 

60 

127.9 

34. 1 

420.7 

24-81-94 

105 

71.0 

41.9 

435.5 

-15.0 

24-81-94 

150 

50.0 

35.2 

436.5 

24-81-94 

side- on 

75.5 

40.2 

587.2 

24-81-94 

stag 

-/93 . 5 

38.5 

676.4 

24-81-95 

30 

172.2 

38.8 

602.6 

24-81-95 

75 

96.6 

32.7 

314.6 

24-81-95 

120 

59. 0 

40.7 

455.7 

-5.0 

24-81-95 

165 

35.5 

34.7 

454.4 

24-81-95 

side- on 

76.0 

40.8 

589.8 

24-81-95 

stag 

180.0/93.5 

38.9 

677.2 

24-81-96 

45 

156.7 

40.0 

499.2 

24-81-96 

90 

88.4 

33.2 

384.1 

24-81-96 

135 

47.0 

40.3 

467.3 

24-81-96 

180 

85.0 

33.6 

481.0 

24-81-96 

side-on 

76.5 

39.6 

587.1 

24-81-96 

stag 

177.5/93.0 

38.7 

677.8 

24-81-97 

0 

297.6 

46.7 

1,119.4 

24-81-97 

45 

219.4 

35.2 

760.0 

24-81-97 

90 

121.8 

39.0 

447.3 

-23.0 

24-81-97 

135 

73.0 

38.8 

524.6 

-25.0 

24-81-97 

side-on 

112.0 

48.5 

905. 1 

24-81-97 

stag 

287.6/146.5 

47.7 

1,112.4 

24-81-98 

15 

295.1 

49.7 

1,109.4 

24-81-98 

60 

193.2 

35.0 

567.0 

24-81-98 

105 

104.5 

51.6 

585.1 

-35.0 

24-81-98 

150 

60.0 

38.7 

558.4 

24-81-98 

side-on 

114.5 

47.5 

919.3 

24-81-98 

stag 

295.9/152.0 

46.8 

1,127.3 

24-81-99 

30 

265.1 

49.7 

952.0 

24-81-99 

75 

155.3 

35.4 

409.4 

-38.5 

24-81-99 

120 

84.5 

53.0 

570.7 

24-81-99 

165 

50.0 

38.8 

603.9 

24-81-99 

side-on 

112.0 

48.1 

902.4 

24-81-99 

stag 

288.2/146.5 

41.5 

1,102.4 

24-81-100 

45 

240.6 

44.5 

718.8 

24-81-100 

90 

124.2 

35.4 

453.4 

24-81-100 

135 

65.0 

51.7 

524.1 

-40.0 

24-81-100 

180 

124.5 

37.4 

620.0 

24-81-100 

side-on 

110.3 

50.8 

903.1 

24-81-100 

stag 

281.1/147.0 

44.4 

1,104.5 

*Positive  impulse  given  for  10  milliseconds. 
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TEST:  2-0  CYLINDER 


TEST:  2-0  CYLINDER 


TEST:  2-0  CYLINDER 


Figure  6,  Pressure-time  records  for  0  degrees. 
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PRESSURE,  kPa  PRESSURE 


TEST:  2-D  CYLINDER 


TIME,  MSEC 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  7. 


Pressure- time  records  for  15  degrees. 
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PRESSURE,  kPa  PRESSURE 


TEST:  2-0  CYLINDER 


TIME,  MSEC 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  8.  Pressure-time  records  for  30  degrees. 
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B.  Results  for  Stations  45-75 


Stations  for  this  area  of  the  model  show  initial  peak  pressures  less  than 
full  reflected  pressures.  The  peak  value  becomes  less  with  increasing 
station  (angle).  At  Station  75  the  initial  pressure  is  65-75  percent  of  the 
initial  pressure  at  Station  45. 

The  general  shape  of  the  records  are  very  similar  for  all  these  stations. 
At  Station  75  the  peak  is  followed  by  rarefactions  that  began  to  dip  toward 
zero  overpressure.  This  is  seen  at  the  higher  two  input  pressure  levels. 

C.  Results  for  Stations  90-120 


At  Station  90,  for  the  75.9  kPa  input  level  pressure,  a  negative  pressure 
occurs  at  about  1  ms.  The  rarefaction  decay  plus  the  vortex  action  causes 
this  pressure  decrease.  This  has  not  occurred,  however,  at  the  lowest  input 
pressure  level.  The  values  for  these  negative  pressures  are  listed  in  Table 

2. 


The  initial  pressure  peak  has  continued  to  decrease  in  amplitude  until 
at  Station  120  the  initial  peak  value  is  no  longer  the  largest  value. 

D.  Results  for  Stations  135-180 


The  negative  pressure  is  absent  on  records  from  Station  135  at  the  two 
lower  input  pressure  levels.  At  Station  150  the  negative  pressure  is  absent 
at  the  highest  input  pressure  also.  The  remaining  stations  show  records  that 
look  similar  to  the  side-on  input  records.  The  rounded  front  records,  Figure 
17,  show  a  two  part  rise  to  the  maximum  value  (neglecting  the  small  reflection 
at  about  1  ms).  The  two  parts  are  caused  by  different  arrival  times  of  the 
diffracted  waves  around  the  cylinder. 

For  example,  the  first  rise  at  the  165°  position  is  the  initial  shock 
front  passing  the  transducer  position.  The  second  rise  in  pressure  is  from 
the  interaction  of  the  two  parts  of  the  initial  shock  meeting  at  180°  and 
passing  back  upstream  over  the  165°  position.  At  Station  180  the  arrival 
times  for  the  two  parts  of  the  initial  shock  (around  each  side  of  the  cylinder) 
are  equal  and  a  single  shock  reflection,  building  up  to  the  rounded  maximum, 
is  the  result. 

Tables  2  and  3  summarize  the  test  results. 


IV.  ANALYSIS 

The  analysis  will  treat  three  topics.  The  first  section  will  show  a  method 
for  determining  the  free-air,  blast  equivalent  of  the  shock  wave  produced 
by  the  shock  tube  for  the  experimental  loading.  The  second  section  will 
describe  a  way  to  calculate  the  coefficient  of  drag  as  a  function  of  time  for 
the  blast  loads  measured  on  the  cylinder.  The  third  section  will  be  a  presen¬ 
tation  of  pressure  coefficients  (P  /P  )  versus  angle  of  incidence  for  the 
three  input  pressures. 
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TABLE  3.  INPUT  SHOCK  WAVE  PARAMETERS 
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PRESSURE,  kPo  PRESSIRE,  kPo 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  9.  Pressure -time  records  for  45  degrees. 
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PRESSURE,  kPo  PRE55WE,  kPfl 


TEST:  2-0  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  10.  Pressure-time  records  for  60  degrees. 
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kPa  PRESSURE,  kP* 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  11.  Pressure-time  records  for  75  degrees. 
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TEST:  2-0  CYLDfcER 


TEST:  2-0  CYLINDER 


TEST:  2-0  CYLINDER 


Figure  12.  Pressure-time  records  for  90  degrees. 
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PRESSURE,  kPa  PRESSURE,  kPo 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  13.  Pressure-time  records  for  105  degrees. 
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PRESSURE,  kPo  PRESSURE,  kPa  PRESSURE 


168 


TEST:  2-0  CYLINDER 
SHOT:  24-81-81 
STATION:  128 


TEST:  2-0  CYLINDER 


TEST:  2-0  CYLINDER 


Figure  14,  Pressure-time  records  for  120  degrees. 
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PRESSURE,  kPa  PRESSURE 


TEST:  2-0  CYLINDER 


TIME,  MSEC 


TEST:  2-0  CYLINDER 


TEST:  2-D  CYLINOER 


Figure  15.  Pressure-time  records  for  135  degrees. 
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PRESSURE,  kPo  PRESSURE,  kPa  PRESSURE,  kPa 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST;  2-D  CYLDOER 


Figure  16.  Pressure-time  records  for  150  degrees. 
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PRESSURE,  kPo  PRESSURE,  kPa 


TEST:  2-0  CYLIM>ER 


TEST:  2-D  CYLINDER 


TEST:  2-0  CYLINDER 


Figure  17.  Pressure-time  records  for  165  degrees. 
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TEST:  2-0  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  18.  Pressure-time  records  for  180  degrees. 
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A.  Equivalent  Charge  Weight 


Cube  root  scaling  allows  the  blast  parameters  from  one  high  explosive 
yield  of  TNT  to  be  found  for  another  yield.  For  the  same  atmospheric  test 
conditions,  the  scaling  relationships  are  given  by  Equation  1  where  the 
scaling  is  from  charge  (1)  to  charge  (2). 


(1) 


D,  TA,  t,  I,  and  W  are  the  station  distance,  time  of  arrival,  positive 
duration,  positive  impulse  and  charge  mass  of  the  explosive.  Subscript  (1) 
parameter  values  are  taken  from  Reference  12  and  are  listed  in  Table  4  for 
Cases  I,  II,  and  III. 


The  equivalent  charge  mass  to  be  found  may  be  obtained  by  rewriting  a 
portion  of  Equation  1  as  Equation  2. 


W 


2 


where 


(2) 


W2  is  the  equivalent  weight  of  TNT  needed  to  reproduce  a  blast  wave  with  the 
same  side-on  overpressure  and  impulse,  I?.  Table  4  lists  these  average 
values  for  each  set  of  test  shots  for  the  shock  tube  cases.  The  values  used 
for  I,  are  listed  for  Cases  I,  II,  and  III  correspond  to  average  side-on 
pressure  values  obtained  during  the  shock  tube  tests. 

After  equivalent  values  of  W0  are  calculated  (last  column  of 
Table  4) the  remaining  parameter  of  distance,  arrival  time,  and  positive 
duration  may  be  calculated  by  use  of  Equation  1  above.  For  example,  for 
free-air,  a  blast  equal  in  pressure  (75.88  kPa)  and  positive  impulse  (1118.6 
kPa-ms)  to  the  middle  group  of  shock  tube  shots  (24-81-93  to  24-81-96)  would 
be  produced  by  an  equivalent  mass  of  4,244  kg  of  TNT.  The  pressure  would 
occur  at  a  distance  of  48.2  m  from  the  charge  center  of  detonation.  It 
would  arrive  63.2  ms  after  detonation  with  a  positive  duration  of  34.8  ms  and 
with  the  required  1118.6  kPa-ms  of  positive  impulse. 
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Samuel  Gladstone  and  Philip  J.  Dolan- Editor s,  " The  Effects  of  Nuclear 
Weapons , "  Dept,  of  Army  Pamphlet  No.  50-3-  Hdq.  Dept,  of  Army,  March  1977 . 

12 

"Structures  to  Resist  the  Effects  of  Accidental  Explosions, "  TM  5-1300 , 
Dept,  of  Army,  June  1969. 
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TABLE  4.  FREE-AIR  BLAST  PARAMETERS  FOR  TNT  EQUIVALENT 
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Table  4  above  summarizes  the  calculations  for  the  three  blast  overpressure 
levels  used  during  the  cylinder  tests. 

B.  Coefficient  of  Drag 

1 3 

It  is  customary  to  present  loading  data  for  a  test  object  in  a  form 
such  that  a  coefficient  of  drag,  Cp,  might  be  found  for  the  object.  It  will 
be  assumed  that  the  net  horizontal  load  (in  the  diffraction  as  well  as  the 
drag  phase)  across  the  test  structure  can  be  found  from  the  appropriate 
measured  pressure-time  profile  multiplied  by  a  normal  projected  area.  The 
sum  of  these  loads  over  the  normal  area  will  be  the  total  horizontal  load 
exerted  on  the  object  by  the  pressure  from  the  blast  wave.  Equation  3 

expresses  these  relationships  for  the  coefficient  of  drag,  C  . 

o  D 

180 

26A^^  p(0) 

CD  qA  qA  5  ^ 

where  the  drag  force,  F,  is  obtained  from  the  summation  of  the  net  pressure 
difference  across  the  cylinder  for  the  projected  normal  surface,  A.  6A  (See 
Table  5.)  is  the  incremental  projected  normal  surface  calculated  for  7.5°  each 
side  of  a  transducer  position.  P(0)  is  the  pressure  as  a  function  of  angle. 

Cp  changes  as  a  function  of  time  with  pressure  changes  in  the  dynamic  pressure, 
corresponding  to  variations  in  the  free-field  blast  wave. 

Equation  4  gives  the  relationship  used  to  calculate  q. 


q  = 


2.5  P 


7P  +  P 
1  s 


(4) 


where  P^  is  the  side-on  overpressure  and  P  is  the  ambient  pressure.  Equation 
4  is  strictly  correct  only  at  the  front  of  the  blast  wave  but  will  be  used 
throughout  the  entire  blast  wave's  positive  duration  for  calculations  of  q 
as  a  function  of  time.  This  method  agrees  well  with  the  exponential  decay  of 
q  as  given  in  Reference  14.  The  method  used  here  was  considerably  more 
convenient  to  the  data  processing  method  used  than  the  more  exact  method 
given  in  Reference  14. 
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Sighard  F.  Hoerner 3  Drag^  Published  by  author 3  148  Busteed 

Drive 3  Midland  Park 3  NJ3  1965. 


"Design  of  Structures  to  Resist  Effects  of  Atomic  Weapons 3  Weapons  Effect 
Data3  "  Dept,  of  Army  Technical  Manual  TM  5-856-1 3  November  1960. 
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Figure  19  shows  the  force-time  curves  computed  from  the  side-on  blast 
wave  after  being  multiplied  by  the  total  normal  area/length,  m.  The  pressure 
records  used  were  smoothed  through  the  small  pressure  peaks  propagated 
upstream  from  the  cylinder.  This  was  done  to  better  represent  undisturbed 
free-stream  flow. 

Table  5  lists  the  pertinent  parameters  needed  for  the  calculation  of  total 
drag  force,  F,  as  a  function  of  time.  This  function  is  displayed  in  Figure 
20.  Substituting  the  data  from  Figure  19,  force  versus  time,  (from  q  versus 
time)  and  Figure  20  (for  drag  F,  versus  time)  into  Equation  3,  the  coefficients 
of  drag  versus  time  can  be  computed  for  the  cylinder.  This  was  done  and  is 
presented  in  Figure  21. 

It  should  be  noted  from  Figure  21  that  the  initial  loading  of  the  cylinder 
during  the  diffraction  phase  caused  coefficients  of  drag  several  times  the 
average  value  after  the  diffraction  phase.  A  second  point  of  interest  is  that 
Figures  21-A  and  21-B  show  negative  values  for  the  coefficients  for  the 
interval  1  to  1.5  ms.  The  remainder  of  the  records  show  oscillations  in  the 
drag  coefficients. 

Vortex  shedding  (Reference  15,  16,  and  17)  could  account  for  the 
oscillations  present.  A  meaningful  average  value  is  a  bit  difficult  to 
determine;  however,  Table  6  and  Figure  22  attempt  to  relate  an  average  early 
time  coefficient  (2-6  ms)  for  the  present  work  and  that  reported  in  References 
15  and  16  with  steady  state  values  from  Hoerner  (Reference  13).  Table  6  lists 
the  average  coefficients  determined  for  early  times  after  the  diffraction 
phase  as  a  function  of  Reynolds  number  (using  diameter  of  the  cylinder)  and 
flow  Mach  number  behind  the  shock  front. 

The  present  results  and  the  shock  tube  results  reported  in  Reference  15 
tend  to  cluster  about  the  upper  curve  of  Figure  22  indicating  Reynolds  numbers 
below  transition  flow.  Whereas,  the  field  data  reported  from  the  Dice  Throw 
Event  (Reference  16)  clusters  about  the  lower  curve  of  Figure  22,  indicating 
flows  above  transition  Reynolds  numbers.  To  predict  accurately  the  coefficient 
of  drag,  it  is  necessary  for  the  test  structure  to  be  clearly  in  one  or  the 
other  of  the  two  flow  regions.  In  the  Mach  number  region  0.45  to  0.50,  one’s 


15 

Valerie  C.  Martin,  K.F.  Mead ,  and  J.E.  Uppard ,  "The 
Cylinder  in  a  Shook  Wave,"  AWRE  Report  No.  0-34/67, 


Drag  on  a  Circular 
May  1967. 


A.W.M.  Gibb  and  D.A.  Rill,  "Free  Flight  Measurement  of  the  Drag  Forces  on 
Cylinders  in  Event  Dice  Throw, "  Suffield  Tech  Paper  No.  453,  February  1979. 
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FORCE,  kN/m  FORCE,  kN/m  FORCE,  kN/m 


TEST:  2-D  CYLINDER 
SHOT:  24-8I-C88  TO  92) 

FORCE  =  qA 


TEST:  2-D  CYLINDER 


B.  Input  of  75.9  kPa 

TEST:  2-D  CYLINDER 


C.  Input  of  112.2  kPa 

Figure  19.  Net  horizontal  force-time  functions 
computed  from  the  side-on  pressure¬ 
time  records. 
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TABLE  5.  DRAG  PARAMETERS 


Ambient 

Positions  Pressure 

kPa 


0  - 

180a 

102.7 

15  - 

165 

102.7 

30  - 

150 

102.7 

45  - 

135 

102.7 

60  - 

120 

102.7 

75  - 

105 

102.7 

P  = 
s 

42.25 

kPa, 

Pj  =  102.7 

kPa, 

T1 

0  - 

180 

102.9 

15  - 

165 

102.9 

30  - 

150 

192.9 

45  - 

135 

102.9 

60  - 

120 

102.9 

75  - 

105 

102.9 

P  = 
s 

75.88 

kPa, 

Pj  =  102.9 

kPa, 

T1 

0  - 

180 

102.9 

15  - 

165 

102.9 

30  - 

150 

102.9 

45  - 

135 

102.9 

60  - 

120 

102.9 

75  - 

105 

102.9 

P  = 

112.2 

kPa, 

P.  =  102.9 

kPa, 

T\ 

Area/Length  Shots 

m 


. 0132b 

24-81-88  - 

24-81-92C 

.0256 

24-81-90  - 

24-81-91 

.0230 

24-81-91  - 

24-81-90 

.0187 

24-81-88  - 

24-81-88 

.0132 

24-81-90  - 

24-81-91 

.0077 

20. 9°C 

24-81-91  - 

24-81-90 

.0132 

24-81-93  - 

24-81-96 

.0256 

24-81-94  - 

24-81-95 

.0230 

24-81-95  - 

24-81-94 

.0187 

24-81-93  - 

24-81-96 

.0132 

24-81-94  - 

24-81-95 

.0077 

24-81-95  - 

24-81-94 

21 . 0°C 

.0132 

24-81-97  - 

24-81-100 

.0256 

24-81-98  - 

24-81-99 

.0230 

24-81-99  - 

24-81-98 

.0187 

24-81-97  - 

24-81-97 

.0132 

24-81-98  - 

24-81-99 

.0077 

21 . 2°C 

24-81-99  - 

24-81-98 

cL 

Area  includes  factor  of  two  to  include  both  instrumented  and 
uninstrumented  sections  of  model . 


Normal  projected  areas  were  calculated  for  7.5  degrees  each  side  of 
transducer  position  except  for  75  and  105  degree  positions.  Area  to 
the  90  degree  point  was  included  for  these  two  positions. 

c  2 

Total  projected  normal  area,  A,  for  all  shots  is  0.1016  m  /length. 
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FORCE,  kN/m  FORCE,  kN/m  FORCE,  kN/m 


TEST:  2-D  CYLINDER 


TIME,  MSEC 

A.  Input  of  42.3  kPa 


TEST--  2-D  CYLINDER 


TIME,  MSEC 

B.  Input  of  75.9  kPa 


TEST:  2-D  CYLINDER 


C.  Input  of  112.2  kPa 

Figure  20.  Net  horizontal  force/length  of 
model  as  a  function  of  time. 
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COEF  DRAG  COEF  DRAG  COEF  DRAG 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  21.  Coefficients  of  drag  computed 
for  the  net  horizontal  loads. 
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TABLE  6.  COEFFICIENT  OF  DRAG,  2-D  CYLINDER  PARAMETERS 
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^Average  values  of  Reynolds  No. ,  Mach  No. ,  and  coefficient  of  drag  over  2-6  millisecond  of 
record  time. 
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Figure  22.  Average  coefficient  of  drag  as  a  function  of  Mach  number 


predictions  might  range  from  the  low  curve  value  for  Cp  of  0.34  to  1.20  for  the 
high  curve.  Great  care  is  needed  for  accurate  prediction  of  drag  coefficients 
for  test  structures. 

C.  Pressure  Coefficient  versus  Angle 

Also  of  interest  in  determining  the  response  of  a  structure  is  knowing  the 
enhancement  of  the  side-on  pressure  as  a  function  of  angle  of  incidence.  The 
reflected  pressures  P  divided  by  side-on  pressures  P  taken  from  Table  2  are 
plotted  in  Figure  23.  Values  of  P  /P  overlap  at  45  degrees,  and  beyond  45 
degrees  can  be  represented  by  a  single  curve. 


V.  SUMMARY  AND  CONCLUSIONS 

A  two-dimensional  non-responding  cylinder  was  tested  in  the  BRL  57.5  cm 
shock  tube  at  nominal  overpressure  levels  of  42.3,  75.9,  and  112.2  kPa. 
Pressure-time  loading  records  were  obtained  at  15°  increments  of  the  cylinder 
for  each  of  the  pressure  levels  for  a  decaying  input  shock  wave.  These  blast 
parameters  scaled  to  those  that  would  be  produced  by  detonation  of  712, 

4244,  and  11898  kg  of  TNT  high  explosive. 

Total  loads  from  the  pressure  loading  were  obtained  by  summation  of  the 
pressure  over  the  projected  normal  area  of  the  cylinder.  Coefficients  of 
drag  as  a  function  of  time  were  computed  from  the  net  horizontal  drag  forces 
and  the  free-field  input  waves  from  which  the  dynamic  pressures,  q,  had  been 
calculated. 

The  resulting  coefficients  of  drag  curves  were  presented  for  the  three 
pressures  tested.  A  great  variation  of  several  times  for  the  coefficient 
was  seen  from  the  values  of  the  diffraction  portion  to  the  values  of  the 
semi-steady  state  portion  of  the  curves.  An  early  time  (2-6  ms)  average 
coefficient  was  compared  to  the  steady  state  coefficients  listed  in  Hoerner 
(Reference  13). 

It  is  seen  from  these  comparisons  that  the  flow  region  to  which  the 
cylinder  is  exposed,  has  to  be  clearly  defined  in  order  to  make  accurate 
predictions  of  transient  drag  coefficients. 
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APPENDIX  A 

DRAWINGS  FOR  CYLINDER 
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NOTES  : 

1.  STEEL  PIPE  TO  BE 
SUPPLIED. 

2.  CUT  TO  LENGTH, 
DRILL  &  TAP 
HOLES. 


DRILL  &  TAP  4  EA.  3/8-16  NC, 
1"  -  DEEP  ON  3"  B.C. 

90°  APART 


DRILL  &  TAP 
1/2-20  NF-2A  FIT 
SPACED  AT  45° 


Figure  A-l.  Sketch  of  two-dimensional  model  cylinder. 
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NOTES  : 

1.  MTL  -  STEEL 

2.  MACHINE  CENTER 
POINT  ON  PLUG. 

3.  MAKE  ONE. 


Figure  A-2.  Bottom  plug  for  model. 
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NOTES: 

1.  MOUNTING  PLATE  TO  BE  FURNISHED. 

2.  TRY  BOTTOM  PLUG  FOR  FIT. 

3.  FURNISH  3/8-2  3/4  LONG  MOUNTING  BOLTS. 

4.  SCRIBE  MOUNTING  PLATE  AT  FRONT  -  0°. 


Figure  A-3.  Assembly  drawing. 
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APPENDIX  B 

PRESSURE-TIME  RECORDS 
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PRESSURE,  kPo  PRESSURE,  kf>a 


TEST:  2-D  CYLINDER 
SHOT;  24-81-87 
STATION:  0 


TIME,  MSEC 


TEST:  2-D  CYUM)ER 


TEST:  2-D  CYLINDER 


Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa. 
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PRESSURE,  kPa  PRESSURE,  kPa  PRESSURE, 


168 


TEST:  2-0  CYUM>ER 
SHOT:  24-81-87 
STATION:  135 


O 

Sc 


.  60 


TEST:  2-0  CYLINOER 


TEST:  2-0  CYLINDER 


Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa*  (cont'd) 
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PRESSURE,  kPa  PRESSURE 


TEST:  2-0  CYUMJER 
SHOT:  24-81-88 
STATION:  8 


TIME,  MSEC 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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PRESSURE,  kPo  PRESSURE,  kPo  PRESSURE,  kPo 


TEST:  2-D  CYLINDER 


100 


TEST:  2-D  CYLINDER 
SHOT:  24-61-69 
STATION:  105 


80 


60 

40 


20 


2 


4 


6 


1 


8 


J. 

10 


TIME,  MSEC 


TEST:  2-D  CYLINDER 


Figure  B-l.  Pressure- time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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pressure: 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


TEST:  2-D  CYLINDER 


Figure  B- 1 .  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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PRESSURE,  kPo  PRESSURE, 


100 


TEST:  2-D  CYLINDER 
SHOT:  24-81-89 
STATION:  SIDE-ON 


80 

60 


TEST:  2-D  CYLINDER 


Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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Figure  B-l.  Pressure-time  records  from  cylinder  for 

input  overpressure  of  42.3  kPa.  (cont'd) 
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Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  fcont'd) 
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Figure  B- 1 .  Pressure-time  records  from  cylinder  for 

input  overpressure  of  42.3  kPa.  (cont'd) 
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Figure  B- 1 .  Pressure-t ime  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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Figure  B-l.  Pressure-t ime  records  from  cylinder  for 

input  overpressure  of  42.3  kPa.  (cont’d) 
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Figure  B-l.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  42.3  kPa.  (cont'd) 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa. 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-2.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-2,  Pressure-time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont’d) 
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Figure  B-2.  Pressure- time  records  from  cylinder  for 
input  overpressure  of  75.9  kPa.  (cont'd) 
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Figure  B-3.  Pressure- time  records  from  cylinder  for 
input  overpressure  of  112.2  kPa. 
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Figure  B-3.  Pressure-time  records  from  cylinder  for 

input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-t ime  records  from  cylinder  for 

input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-time  records  from  cylinder  for 

input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-time  records  from  cylinder  for 
input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-t ime  records  from  cylinder  for 
input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-t ime  records  from  cylinder  for 

input  overpressure  of  112.2  kPa.  (cont'd) 
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Figure  B-3.  Pressure-time  records  from  cylinder  for 

input  overpressure  of  112.2  kPa.  (cont'd) 
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5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) _ 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


